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ABSTRACT 

A computer model simulating the seasonal variations of 
mixed layer nutrient concentrations, phytoplankton biomass 
carbon, and herbivorous zooplankton biomass carbon was devel- 
oped, The simulation was generated using an annual cycle of 
four environmental parameters! (1) Incident solar radiation, 
(2) upwelling velocity, (3) mixed layer depth, and (4) mixed 
layer temperature. Simulation results were compared with 
nutrient and zooplankton biomass data collected on a series 
of seven cruises made in central Monterey Bay from February 
through December, 1974. Both observed and simulated zoo- 
plankton stocks were characterized by two distinct maxima. 

2 

The initial peak (1.05 gC/m ) occurred in late July and was 

followed by a decline in populations through the month of 

August. During the fall and early winter, zooplankton biomass 

2 

increased rapidly to an overall maximum of 1.85 gC/m . 
Individual environmental parameters were tested to ascertain 
their importance in controlling simulation results. Phyto- 
plankton stocks were influenced principally by changes in 
incident radiation, whereas temperature variations produced 
the most significant fluctuations in zooplankton biomass. 
Simulation responses suggest that upwelling, in addition to 
providing nutrients for primary production, enhances zoo- 
plankton productivity by bringing colder deep water to the 
surface, thereby reducing zooplankton respiration requirements. 
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I, INTRODUCTION 



A . PURPOSE 

There presently exists a large body of detailed informa- 
tion about the biological, physical, and chemical processes 
which occur in the marine environment. These processes act 
in concert to control the plankton populations of the sea. 
The analysis and interpretation of these environmental and 
biological factors has been significantly enhanced by the 
development of computer simulation techniques. The aim of 
computer modeling is to reproduce the dynamic response of a 
biological ecosystem to the complex mix of interacting 
variables which act on and within that system. 

The purpose of this research was threefold: 

(1) to formulate a model of a planktonic ecosystem in a 
region of upwelling, 

(2) to develop a computer simulation of the model which 
generates seasonal changes in mixed layer nutrient 
concentrations, phytoplankton biomass carbon, and 
zooplankton biomass carbon, and 

(3) to test the accuracy of the simulation results by 
comparison with observed data. 

B. CONCEPTUAL MODEL 

The first step in the modeling process consists of 
reducing the natural system being simulated to a simplified 
conceptual model. From this conceptual model, processes and 
rates are formulated into mathematical equations. The model 
is completed by translating the equations into the proper 
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format. The conceptual model considered here describes the 
relationships between factors controlling biological 
productivity (see Fig. 1). 

1 , Production 

The synthesis of organic compounds from the inorganic 
constituents of the sea by the activity of living organisms 
is termed ’’production" (Tait and DeSanto, 1972), This 
synthesis is carried out almost exclusively by the photo- 
synthetic activity of marine plants. The raw materials for 
photosynthesis are water, carbon dioxide and various nutrients, 
mainly inorganic ions of nitrogen CNO^, NO 2 , NH^) and 
phosphorus (P0“, HPO^, H2P0~) . The chlorophyll-containing 
plants utilize light energy (nhy) to synthesize higher energy 
organic compounds from the lower energy inorganics as follows; 

6 nhy + 6 CO 2 + 6 H 2 O 6 CH 2 O + QO^ (1) 

This synthesis is referred to as "gross primary production." 

A certain percentage of this synthesized material is 
broken down by the plant cells to provide energy to carry on 
physiological processes. This loss is termed respiration. 

The remainder of the synthesized material is available for 
cell growth and replication and is referred to as "net 
primary production." It is the net primary production of 
phytoplankton that is available for consumption by the 
herbivorous zooplankton, 

Grazing by herbivores on the phytoplankton leads to 
"secondary production." Herbivores in turn are a food source 
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for first level carnivores ("tertiary production"), who them- 
selves become prey for other carnivores, and so on. This 
type of organization of a biological community is referred 
to as a "trophic structure" or "food chain," At each level 
losses of biomass may occur from respiration, excretion, 
predation, or natural mortality. In reality the relation- 
ships are more complicated than this simple "food chain" 
concept indicates. Since one organism may take food from 
more than one trophic level, a more correct structure would 
define an intricate "food web," 

In spite of its over-simplification, the "food chain" 
model depicts certain important characteristics of the marine 
planktonic ecosystem. First, the path of energy in the form 
of organic material is part of a closed cycle in which 
nutrients are assimilated, transferred, and regenerated; and, ~ 
since water and carbon dioxide are abundant in the sea, the 
factors which potentially limit photosynthesis are the 
availability of nutrients and light energy to the phytoplankton 
(Ryther, 1963), and the temperature of the medium (Parsons 
and Takahashi , 1973). 

2 . Light 

The light energy available to the phytoplankton is 

dependent on two parameters, the amount of radiation incident 

2 

on the sea surface (usually expressed as g-cal/cm /min) and 
the transparency of the water (expressed as the extinction 
coefficient, k). Ryther (1956) observed that photosynthesis 
increases linearly with light intensity up to a saturation 
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point. Beyond this level greater illumination produces no 
further increase in photosynthesis and eventually results in 
an inhibition of the photosynthetic rate. 

Light penetrating into the water is absorbed and 
scattered selectively. The intensity is attenuated exponen- 
tially with depth and can be expressed in terms of the 
extinction coefficient (k): 




( 2 ) 



= radient energy at the surface 
I = radient energy at depth z 

2j 

k = light extinction per meter depth. 

The transparency of the water is dependent on the 
concentration of phytoplankton and the presence of particulate 
or dissolved substances which absorb light. If we disregard 
selective absorption of light of different wavelengths and 
consider the total visible spectrum, we may say that for 
clear ocean water, k equals 0,04 - 0.05 (Ryther, 1963). Riley 
(1956) determined empirically the relationship between 
phytoplankton density and transparency. He assumed that 
absorption by material other than water and plants was 
negligible. His relationship is: 

k = 0.04 + 0.0088 Cl + 0,054 Cl^^^ (3) 



where Cl is chlorophyll concentration (yg/1) and 0.04 is 

the extinction coefficient for pure sea water. It is apparent 

from this formula that one mechanism for limiting phtoplankton 
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productivity is self r-shading . Higher concentrations of algae 
produce a higher extinction coefficient and, from equation (2) 
reduce the light intensity in the water column. The reduced 
intensity lowers the rate of photosynthesis. 

Measurements of the photosynthetic rate under 
conditions of bright daylight indicate that radiation intensity 
near the surface is at or above the saturation point for most 
phytoplankton. Under these conditions maximum productivity 
occurs at a depth, usually between 5 20 m. Below the depth 

of maximum productivity the photosynthetic rate decreases as 
the light intensity is attenuated. A point is reached where 
the rate of photosynthesis is just sufficient to balance the 
rate of breakdown of organic materials by respiration. This 
point is referred to as the ’’compensation depth" (Tait and 
DeSanto, 1972). 

A second level below the compensation depth is the 
"critical depth." The concept of the "critical depth" was 
first suggested by Grann and Braarud (1935) and developed into 
a mathematical model by Sverdrup (1953). The critical depth 
is the level at which the total respiration of the phyto- 
plankton in the water column just equals the total photo- 
synthesis. The distance between the compensation depth and 
the critical depth depends on the proportion of the 
phytoplankton stock above and below the compensation depth. 

If the standing stock of phytoplankton is to increase, 
its total phtosynthesis must be greater than its total 
respiration. This is possible as long as the depth to which 
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the plankton Is mixed does not exceed the critical depth. 

When the depth of the mixed layer extends below the critical 
depth (assuming uniform plankton mixing), total losses exceed 
total gains and the stock declines. 

3 . Nutrients 

Another parameter which can critically limit the 
growth of phytoplankton is the concentration of inorganic 
nutrients in the upper, "euphotic" layers of the ocean. Many 
elements such as manganese, silicon, copper, molybdenum, and 
cobalt are essential for the growth of particular species. 
However, ions of nitrogen and phosphorous may be present in 
sufficiently low concentrations that their availability exerts 
a dominant control over production (Tait and DeSanto, 1972). 

The absorption of nutrients in the surface layer by 
phytoplankton can reduce the concentration of these ions to 
levels where further uptake is inhibited. The kinetics of 
this inhibition was investigated by Ketchum (1939). He 
demonstrated that the uptake of nitrate by the diatom 
Nitzschia closterium was concentration dependent over the 
approximate range of 1 - 7 ygr-atom N/1, He also showed that 
phosphate followed a similar inhibition. Dugdale (1967) 
described the rate of nutrient uptaking using Michaelis- 
Menten enzyme kinetics in which: 



^ ^max K + N 
? s 

V = rate of nutrient uptake 

V = maximum rate of nutrient uptake 
max 



( 4 ) 
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K = half saturation constant, i,e. the nutrient 

concentration at which v = V /2. 

max' 

N = nutrient concentration 

Studies conducted by Eppley and Thomas (1969) on two species 
of diatoms indicate that the half-saturation constants for 
growth and nutrient uptake are very similar. This means that 
the expression given for velocity of nutrient uptake can also 

be used to describe the growth rate; 

_ N 

^ max K + N 
s 

y = growth rate 

u = maximum growth rate 

^max ^ 

By plotting N/y versus N, a straight line is obtained with 
the intercept on the abscissa at 

a. Temperature 

Ichimura (cited in Parsons and Takahashi , 1973) 
studied the relationships between phosphate concentration 
and the rates of organic synthesis in natural populations of 
phytoplankton in Tokyo Bay. His results indicate a definite 
Michaelis-Menten type response which is highly temperature 
dependent (see Fig. 2). 

b. Excretion 

A certain fraction of the nutrients absorbed by 

phytoplankton are returned to the medium both in the euphotic 

zone and at depth by zooplankton. According to Harvey (1957), 

nearly all phytoplankton is grazed by herbivorous zooplankton, 

A portion of this grazed algae can be voided in an undigested 

or semidigested state. It is unclear, however, how much 

phytoplankton nutrient elimination by zooplankton is returned 

16 



in an available form in the euphotic zone and how much sinks 
to greater depths before being remineralized. Some authors 
(Gardner, 1937; Riley, 1951; and Harris, 1959) consider 
zooplankton to be a significant source of useable nutrient 
in the euphotic zone, while others believe this source to be 
negligible CSteele, 1959; Rigler, 1971). Most recent work 
seems to favor the higher estimates, A distribution of 
dietary phosphorus in Calanus during active feeding (April) 
has been given by Butler et al . (1970) as 17.2% retained for 
growth, 23.0% voided as fecal pellets and 59.8% excreted as 
soluble phosphorus (organic and inorganic). These estimates 
indicate that over 80% of the phosphorus ingested by Calanus 
gets back into the environment and nearly 60% is soluble 
phosphorus. The ability of marine algae to utilize soluble 
organic phosphorus was demonstrated by Kuenzler (1965). 

The importance of zooplankton as a source of 
phytoplankton nutrients was verified in two studies cited by 
Parsons and Takahashi (1973). Cushing (P. 139) "showed that 
during ten weeks of the spring phytoplankton/zooplankton bloom 
in the North Sea, inorganic phosphate did not decrease below 
0.6 yg-atom/1 due to zooplankton excretion. However, Antia 
et al. (1963) observed a decrease in phosphate to 0.1 yg-atom/1 
after two weeks of phytoplankton growth in the absence of 
zooplankton," suggesting the important role of zooplankton as 
a source of phytoplankton nutrients. 

The nitrogen and phosphorous cycles in the sea 
follow basically the same course, However, the remineralization 
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of nitrogen occurs at a much slower rate than phosphorus. 

In the experiment conducted by Antia et al. (1963) half the 
phytoplankton phosphorus was remineralized in 14 days, but 
after 75 days no remineralization of nitrogen was observed. 

A second major difference between the cycles is 
that all three inorganic nitrogen compounds, ammonium (NH^) , 
nitrite (NOg), and nitrate (NO^) can be absorbed by phyto- 
plankton, Absorption is not indiscriminate, however. There 
is a marked preferential uptake of ammonium (Dugdale and 
Goering, 1967). The nitrogen cycle then consists of three 
seperate cycles, one for each of the three nitrogen containing 
ions involved in photosynthesis, and possibly a fourth involv- 
ing urea (Dugdale, personal communication). 

Although a fraction of the nutrients absorbed by 
the phytoplankton during photosynthesis are recycled within 
the euphotic zone, the remainder are regenerated in deeper 
water. The result is a net depletion of nutrients from 
surface waters and an accumulation of nutrients at depth. 

The continuation of plant growth in the sea, therefore, is 
largely dependent upon the various processes which return 
nutrient-rich water from below the euphotic zone, 
c. Mixing and Upwelling 

Turbulent (eddy) diffusion is a primary mechanism 
of injecting nutrients into surface waters of many ocean areas 
Verticil eddies may be generated by a variety of dynamic 
forces including thermal convection, internal waves, friction 
at current boundaries, and wind-mixing. Turbulence may have 
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a variety of effects on biological production. It may 
enhance production by bringing nutrients to the surface, or 
it may reduce production by carrying phytoplankton down 
below the euphotic zone. 

A second major phenomenon operating in more 
localized regions is coastal upwelling. For example, along 
western continental boundaries in the Northern Hemisphere a 
northerly wind stress blowing parallel to the coast causes 
an offshore transport of water with a compensating replacement 
of upwelled waters from 200»'300 m depth (Sverdrup et al , , 
1942). 

An index of the intensity of upwelling has been 
estimated using wind stress data (Wooster and Reid, 1963; 
Bakun, 1973). This *'Ekman transport^' of water away from the 
coast is calculated according to the following equation: 



M = T /f (6) 

= "Ekman transport” normal to the coast 
Tp = Wind stress parallel to the coast 
f = Coriolis parameter 

Using the magnitude of offshore transport it is 
possible to compute the velocity of the associated verticle 
motion. Wooster and Reid (1963) assumed a steady-state 
offshore "Ekman transport” of 10 kg/cm/sec with a coastal 
upwelling band 50 km wide and estimated the compensating 
vertical motion to be about 50 m/month, a resut compatible 
with previous estimates which ranged from 10-20 m/month to 



19 



80 m/month (Smith, 1968). Although the depth from which these 
waters are upwelled is not believed to exceed a few hundred 
meters, the input of nutrients is sufficient to maintain 
the regions of upwelling at a high level of productivity. 

Horizontal advection of surface waters could have 
either a positive or negative effect on nutrient concentration. 
If ^'upstream” or "offshore*’ waters have been depleted below 
local values, advection into the area under consideration 
would reduce the concentration and conversely, higher values 
would raise local concentrations. This process is discussed 
later . 

4 . Grazing 

Another major factor controlling the productivity of 
phytoplankton is grazing by herbivorous zooplankton. Accord- 
ing to a number of authors (Steemann-Nielsen , 1958; Cushing, 
1959; Steele, 1974) grazing by zooplankton is the primary 
factor affecting the dynamics of plankton populations. 

Mullin (1963) identified two methods of describing 
zooplankton grazing activity, by the "filtering rate" and 
the "food intake rate." The "filtering rate" is the rate at 
which water is swept clear of food organisms. The "rate of 
food intake" is a measure of the biomass of food organisms 
ingested per unit time and depends on the filtering rate and 
the concentration and size of the food organisms. Mullin 
(1963) demonstrated that, for four species of copepods, the 
volume swept clear is inversely proportional to food concen- 
tration. Steele (1974) showed that such an inverse relation 
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corresponds to a hyperbolic rate of intake of biomass. As 
the food concentration increases, the rate of biomass 
ingestion approaches a maximum (see Fig, 3), McAllister (1970) 
observed intake rates which approximate the type given by 
Steele, Wash and Bass (1971) expressed this saturation feed- 
ing behavior in terms of a Michaelisr^Menten equation. This 
hyperbolic response suggests two important characteristics of 
herbivore grazing. First, a minimum phytoplankton concentra- 
tion exists at which grazing ceases; and second, as algal 
concentration increases, herbivores asymptotically approach 
a maximum ingestion rate. 

In summary, it is apparent that the productivity of 
a region is dependent on the availability of light and 
nutrients to the phytoplankton, the temperature of the medium, 
and the grazing and excretion of biomass and nutrients by 
the herbivores. The rates at which these processes proceed 
are not simple functions but are dependent on a complex mix 
of physical, chemical, and biological parameters. In order 
to model the dynamics of this regional plankton community, 
therefore, it is essential to utilize a method of computer 
simulation. Simulation provides a technique for integrating 
the many variables affecting productivity into a system of 
interacting components, 
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II . BACKGROUND THEORY 



A . EARLY WORK 

The basic equations of biological growth kinetics were 
developed in the early 1920 -s by A. J, Lotka (1925), 

Volterra (1926) extended Lotka'^s equations to describe predator- 
prey interactions. The Lotkar^Volterra equations (7) consist 
of a pair of differential equations in which the growth of 
each population is limited by two factors, self inhibition 
(rf, r 2 ) and competition (a, b) (Odum, 1971); 



—1 = r N 
df^ ^r^i 



Ki - 



^aN, 



K. 



(7) 



dNo _ « 

dt^ ^2^2 



^2 - ^2 ^ bN^ 

K„ 



^1’ ^2 ~ quantity of species one and two 

a, b = competition coefficient indicating the negative 
effect one species has on the other 

r^^, 1*2 = the intrinsic growth rates of each species 

Kf, ^2 ~ carrying capacity, the largest population of 
each species the environment can support. 

Fleming C1939) used a modified form of the Lotka-- 

Volterra equations to calculate the seasonal change in 

phytoplankton population: 



dt 



P(a 



(b + ct)X 



( 8 ) 
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p 



= Integral phytoplankton population under a unit 
of water surface 



a 



= specific growth rate of plankton 



b 



= initial grazing rate 



Cb+ct) = grazing rate as a function of time 



c 



= rate of grazing rate change 



t 



time. 



His model emphasizes grazing as the factor limiting phyto- 
plankton populations, Patten C1968) pointed out that, 
although the model may be unrealistic in specifying a constant 
growth rate and a linearly increasing grazing rate, the model 
"worked" because it accurately corresponded to observations 
made by Harvey et al, (1935) in the English Channel. 

B. PHYTOPLANKTON MODEL 

Riley (1946) made a significant improvement upon the lack 
of realism in Fleming's model. He defined the change in 
phytoplankton as: 



P = integral phytoplankton under a square unit of 
water surface. 

Pj^ = photosynthetic rate 

R = respiration rate 

G = grazing rate. 

Instead of using a constant photosynthetic rate, Riley 
defined Pj^ as a function of light intensity, nutrient concen- 
tration and verticle turbulence. 




(9) 
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( 10 ) 



?h = (pl^/kz^) (1 ^ (1 ^ N) (1 r. V) 

P = fraction of light used in photosynthesis 

I = surface light intensity 
o 

k = extinction coefficient 
= depth of euphotic zone 
N = nutrient depletion factor 

V = factor expressing the effect of turbulence 
His expression for the respiration rate was defined in terms 
of temperature; 

tT 

R = R e^^ (11) 

o ' ' 

R = respiration rate 

R^ = respiration rate at 0°C 

r = rate of respiration rate change 

T = temperature. 

The grazing rate was assumed to be proportional to the inte- 
gral quantity of zooplankton; 

G = gZ (12) 

G = grazing rate 

g = proportionality constant between grazing rate and 
zooplankton 

Z = integral zooplankton under a square unit of water 
surface . 

C, PHYTOPLANKTON, NUTRIENT MODEL 

Steele (1958) introduced a model based largely on the 
equations developed by Riley, His contribution is significant 
in that it is one of the first "systems" approaches to plankton 
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dynamics, Where Riley treated nutrient concentration as an 
independent parameter, Steele defined nutrients and phyto- 
plankton as two interr-dependent variables. 



To develop a nutrient equation, Steele made certain 
simplifying assumptions about the environment. The ocean 
volume was divided into two thermallyr-stratif ied layers with 
homogeneous mixing in each layer. Mixing across the boundary 
was controlled by an exchange coefficient (m) . 

The system of equations devised by Steele define the rate 
of change of nutrients and phytoplankton in a surface layer 
40m deep; 



m = mixing coefficient governing exchange between the 
upper and lower layers. 

P = phytoplankton concentration in the upper layer. 



Pj^ = photosynthetic rate 
r = phytoplankton respiration rate 

f = grazing rate, a function of zooplankton populations 
z = zooplankton population 

V = ratio of sea water viscosity at 0°C to that at 
ambient temperature 

c = factor for converting carbon to phosphate 
p = nutrient concentration in upper layer 
p^ '= nutrient concentration in lower layer. 




(13) 



= -c(p^ - r)P - m(p , p^) 



(14) 



= phytoplankton concentration in the lower layer. 
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D. REGIONAL STEADY STATE MODEL 

Riley (1965) added to Steele's (1958) model a term 
expressing the regeneration of nutrients voided by zoo- 
plankton. This model was tested against observations made 
during quasi-steady state summer conditions found in various 
regions of the Atlantic and Pacific Oceans. Basic environ- 
mental variables, including phosphate concentration in the 
deep water mass, the rate of mixing between the two layers, 
and water temperature, were varied according to the charac- 
teristics of each area. For each array of environmental 
factors, steady-state concentrations of phytoplankton, 
zooplankton, and phosphate were calculated for the upper 
layer. Riley concluded that the model "provided a reason- 
able fit for most of the observed variations." 

E. SHORT TERM SPATIAL MODEL 

Walsh and Bass (1971) developed a spatial model of an 
upwelling ecosystem which was used for short term analysis. 
Their system of equations related the upwelling of nutrients 
to the nitrogen biomass of phytoplankton, zooplankton, and 
fish. A Michaelis-Menten type expression was incorporated 
into the model to define both the nutrient uptake rate and 
the grazing rate. 

To conclude, the early models of Riley (1946) and Steele 
(1958) established the fundamental principles used in 
simulating plankton dynamics. Riley (1965) demonstrated 
that regional variations in phytoplankton could be 
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realistically predicted for a quasi-steady state summer 
period. Using local environmental conditions as input 
parameters, his model accurately predicted the plankton 
stocks observed in different geographic areas. Walsh and 
Bass (1971) developed a spatial model to analyze a specific 
region of upwelling. Since their model applied to short 
term (two week) simulations, environmental parameters 
were held constant. 
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Ill , THE MODEL 



The simulation model proposed in this thesis defines a 
one dimensional, nonlinear, planktonic ecosystem in a 
centralized region of coastal upwelling (Monterey Bay). 

The simulation output consists of a theoretical annual cycle 
of mixed layer nutrient concentration, phytoplankton and 
zooplankton biomass. Although the model has been developed 
for a localized area it may be considered applicable to a 
broader geographic region, Bolin and Abbott (1962) determined 
that "results obtained Cfor the Monterey Bay area) should 
reflect in a broad way the oceanographic conditions and 
events occurring over an extensive segment of the eastern 
North Pacific." As in Riley's (1965) model, environmental 
conditions control the output of the system. From observed 
changes in the physical environment, the model generates 
values for phosphate concentration and plankton standing 
stocks . 

A , COMPARTMENT MODEL 

The compartment model (Fig, 1) defines the flow of 
nutrients and organic carbon through a three component eco- 
system consisting of inorganic phosphate, phytoplankton, and 
herbivorous zooplankton. In this diagram (method from Patten, 
(1971)) the state variables of the system are denoted by 
blocks, and signal flows between the "compartments" are 
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designated by unidirectional branches. The relations between 
the state variables are expressed as a system of differential 
equations , 

The three equations in Fig, Cl) describe the change in 
phosphate concentration in the mixed layer (XI), in the bio- 
mass of phytoplankton C^2), and in the biomass of herbivorous 
zooplankton (X3), The underlying principle used in the 
derivation of these expressions is the conservation of mass. 
Therefore, the validity of the system of equations depends 
on having all significant transfers of matter accounted for. 

The movement of substance or energy into or through the 
system is termed a "flux” or "flow" and is described as a 
quantity of mass or energy per unit time, e.g., gram calories 
per square centimeter per day or grams carbon per square 
meter per day. In ecological systems the flux between 
compartments is often dependent upon the state variables 
involved. In such cases the flux is defined in terms of a 
"rate of transfer" or "flow rate," A flow rate is a 
quantity transferred per unit quantity of source or receiver, 
per unit time. The distinction between flux and flow rate 
can best be made in terms of units. If a flux is in grams 
carbon per square meter per day, the corresponding rate of 
flow is in units of reciprocal days. 

To differentiate between fluxes and flow rates in the 
model, fluxes have been denoted by descriptive variable names 
(NUT, REGEN, GRAZ, etc,). Flow rates are designated by a 
Greek letter followed by two digits (TAU12, RH030, LAM30 , etc.). 



29 



mTAU*' identifies a transfer between tropic levels, "RHO" 
signifies respiration, and "LAM” represents losses due to 
mortality and carnivorous predation. The first digit follow- 
ing the Greek letter identifies the source compartment and 
the second, the receiver. A zero indicates a flow into or 
out of the system. 

B. ECOSYSTEM INPUTS AND BOUNDARY COxNDITIONS 

Inputs which originate as energy or material sources 
outside the system act as forcing functions. The forcing 
functions which drive the system of equations consist of 
four environmental parameters, incident solar radiation (RADI), 
sea surface temperature (TEMP), a coefficient of water volume 
exchange across the layer (M), and an estimate of upwelling 
velocity OV) . 

Values for incident solar radiation (RADI) were obtained 
from tables developed by Kimball (p, 103, in Sverdrup et. al . , 
1942), The values desired for latitude 36® N. were estimated 
by interpolating between those tabulated for 42° N. and 30° N. 

Temperature data was obtained during a series of seven 
cruises conducted aboard the R/V ACANIA from February through 
November 1974, Estimates of the mixed layer depth (Z), used 
to define boundary conditions, were developed from the 
empirical formula of Robinson and Bauer (1971), This formula, 
currently used by Fleet Numerical Weather Central (FNWC) 
defines the mixed layer depth as "The depth at which the 
temperature of the water is 2° F. (1.11° C) below the surface 
temperature . " 
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Bakun (1973) developed a series of monthly indices 
describing the intensity of large scale, wind induced coastal 
upwelling at selected locations along the west coast of 
North America, The indices are based on calculations of 
offshore Ekman transport derived from monthly mean surface 
atmospheric pressure data. Values of offshore transport for 
1974 (Bakun, Personal Communication) were used to estimate 
upwelling velocity (W) , A maximum velocity of one meter 
per day was assumed to correspond with the maximum index of 
upwelling. A linear interpolation was then applied to other 
monthly indices to arrive at an annual variation of upwelling 
velocity. The coefficient of exchange (M) was assigned a 
constant value of 0,12. Other sources of nutrient input 
(i.e,, river runoff etc.) have been assumed to be negligible. 

As mentioned earlier, the models of Steele (1958) and 
Riley (1965) postulated a two layered ocean with homogeneous 
mixing in each layer. Phytoplankton growth was confined to 
the upper layer and the phosphate concentration was assumed 
to remain constant in the lower layer. Mixing between the 
two layers supplied nutrients to the surface. The present 
model incorporates these assumptions and includes upwelling 
as an additional source of nutrient input. 

C. ECOSYSTEM EQUATIONS 

Having specified the input signals and boundary conditions, 
it is possible to define a system of equations to generate 
the response of the state variables to variations observed in 
the physical environment. 
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1 , Nutrient Equation 



The first equation, expressing the change in phosphate 
concentration in the mixed layer, is given by: 



dXl 

dt 



NUT + REGEN 



UPTAK 



(15) 



The change 
in nutrient 
concentrat-( 
ion with 
time 



The input 
of nutrient 
)from upwell^r 
ing and 
mixing 



The regent- 
eration inr- 
+)put from (-) 
zooplankton 
excretion 



The loss due 
to photosyn- 
thetic uptake 
by 

phytoplankton . 



The expression defining the input of nutrients from the lower 
layer is; 



NUT + (XW/EKL)*(DXlrOQ)) + ( (M/Z ) ♦ (DXl^l ) ) (16) 



NUT 

W 

EKL 

DXl 

XI 

M 



flux of nutrients into the 
surface layer 

upwelling velocity 

depth of Ekman layer 

nutrient concentration in 
the deep layer 

nutrient concentration in 
the mixed layer 

coefficient of exchange 

mixed layer depth 



(yg atom P/l/day) 
Cm/day) 

(m) 

(yg atom P/1) 

(yg atom P/1) 
(m/day) 

(m) , 



The ratio of the coefficient of exchange to the mixed layer 
depth (M/Z) represents the fraction of the upper layer 
transferred to the lower layer in a days time, with a 
corresponding replacement from below. The ratio of the 
upwelling velocity to the depth of the Ekman layer (W/EKL) 
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represents the fraction of the surface water that is trans- 
ported offshore and replaced by upwelled water. 

The second term in the nutrient equation represents 
the input of phosphate from the nutrients voided by 
zooplankton ; 



REGEN = TAU31 * X3 (17) 

REGEN = flux of nutrients from 

zooplankton (,ygt-atom/l/day ) 

TAU31 = rate of input of phosphate 

into the nutrient pool „ 

by zooplankton (yg-atom/l/g C/m^/day) 

X3 = standing stock of 2 

zooplankton (g C/xa). 

The rate at which nutrients, voided by zooplankton, are made 

available for phytoplankton uptake is a function of the 

grazing rate of the zooplankton. 



TAU31 = FAVL * (C/Z) * TAU30 



(18) 



TAU31 = rate of input of phosphate 

to the nutrient pool by 2 

zooplankton (pg-atom/l/g C/m'^/day) 

FAVL = fraction of voided nu- 
trient made available 

for phytoplankton uptake (dimensionless) 

C = converts biomass to 2 

nutrients voided. (yg-atom/l/g C/m ) 



Z = depth of the mixed 

layer 



(m) 



TAU30 = rate at which biomass 

is voided by ^ 

zooplankton (day”"^). 



The rate at which biomass is voided is dependent upon the 



grazing rate (see equation 37). 
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mm 














The final term in the nutrient equation is an 
expression for the loss of phosphate to phytoplankton: 



UPTAK = C * TAU12 ♦ (X2/Z) 



( 19 ) 



UPTAK = flux of nutrients to 
phytoplankton 



(yg atom P/l/day) 



C 



a conversion factor 
relating organic carbon 
to phosphate 



Cyg atom P/l/g C/m^) 



TAU12 = rate of organic pro 



duction by a unit of 
phytoplankton 




X2 = standing stock of 
phytoplankton 



(g C/m^) 



Z = mixed layer depth 



(m). 



2 , Phytoplankton Equation 

The equation for the second state variable, phyto- 
plankton, is an adaptation of Riley's (1946) model. Modi- 
fications have been made on the basis of recent findings, e.g. , 
the effects of a potentially limiting nutrient on photo- 
synthesis have been expressed in terms of Michaelis-Menten 
kinetics (Dugdale 1967); McAllister (1970) demonstrated that 
zooplankton grazing could be defined in terms of saturation 
feeding; and, Walsh and Bass (1971) used a form of the 
Michaelis-Menten expression to describe both the effects of 
nutrient limitation in photosynthesis and the saturation 
response in zooplankton grazing. The approach of Walsh and 
Bass has been followed here. 

The differential equation for the change in phyto- 
plankton biomass carbon is; 
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- PROD 
dt 


rv RESP2 V GRAZ 




(20) 


which in word form 


is ; 






The change 


The pro^ 


The loss 


The loss to 


in phyto- 


duction by 


due to 


grazing by 


plankton (=) 


photosyn- (r^) 


respira- 


(-) zooplankton. 


biomass 


thesis 


tion 




with time 









The expression defining primary production is a 
product of the photosynthetic rate and the phytoplankton 
biomass : 



PROD = TAU12 * X2 



( 21 ) 



PROD = photosynthetic flux to 
phytoplankton biomass 

TAU12 = rate of organic produc- 
tion by a unit of phyto 
plankton 



(g C/m^/day) 



(day”^) 



X2 = standing stock of 
phytoplankton 



(g C/m^). 



Since photosynthesis in the ocean may be limited by temperature, 



light, or nutrients, each of these parameters must appear in 



the rate equation: 



TAU12 = P * RAD * NTLIM (22) 

TAU12 = rate of organic production ^ 

by a unit of phytoplankton (day” ) 

P = photosynthetic constant, 

converts light energy units 
to photosynthetic 2 

potential (g cal/cm /min) 

RAD = mean radiant energy 

available in a water 2 

column of depth (Z) (g cal/cm /min) 

NTLIM = potential photosynthetic 

rate, a function of temperature 
and phosphate concentration (day ‘^). 
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The value for the radiant intensity at depth Z (RADZ) is 
derived from the incident intensity (RADI) as follows: 



(g cal/cm /min) 
2 

(g cal/cm /min) 
(m"^) 



RADZ = RADI * expC-K*Z) (23) 

o 

RADZ = intensity at depth Z 
RADI = incident intensity 
K = extinction coefficient 
Z = mixed layer depth (m) 

To find the mean radiant energy (RAD) between the surface 
and the mixed layer depth (Z), equation eight is integrated 
from zero to Z and devided by Z, 



RAD = (RADI/Z) * 

O 

RAD = (RADI/(r.K*Z)) * (1,0 r- exp(-K*Z)) (24) 

where exp(^<K*Z) = 1,0 at the surface. 



/■ 



exp(^K*Z)dZ 



The extinction coefficient (K) was estimated using Riley's 
(1956) empirical formula; 

K = 0.04 + CO. 0088 * CL) + (0.054 * CL) 2/3 (25) 

K = extinction coefficient (m~^) 

3 

CL = chlorophyll concentration (mg/m ). 

Riley (1959) indicated that carbon : chlorophyll ratios may 
vary from 30 to 128. During periods of algal growth, the 
ratio approaches the lower limit. It has been assumed that 
phytoplankton densities are significant during this rapid 
growth phase. Therefore, the lower ratio of 30 was used. 
Using this ratio, the value lor the chlorophyll concentration 
becomes : 
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CL « 33,33 * C.X2/Z) , 

The function NTLIM CFig, 2) defines the maximum 
photosynthetic rate possible for a given combination of 
temperature and phosphate concentration. The function 
approximates the results of Ichimura (p, 72, in Parson and 
Takahashi, 1971): 



NTLIM = MAXN * (XI ^ X1MIN)/(KN + (XI XIMIN)) (27) 

NTLIM = photosynthetic rate as a 

function of nutrient con^ ^ 

centration and temperature (day" ) 



XI = nutrient concentration 

in the mixed layer (pg atom/1) 

XIMIN =5 minimum nutrient concen^^ 
tration below which 

phytosynthesis will not occur (pg atom/1) 
KN = half saturation constant (pg atom/1). 



The value for MAXN is defined as a function of temperature: 



MAXN = S * TEMP + B (28) 

MAXN = maximum rate of photo^ 
synthesis at a given 

temperature (day" ) 

S = rate of change of MAXN 

with temperature (day /°) 

B = MAXN at 0®C (®C) 

TEMP = mixed layer temperature (®C), 



In summary, the actual photosynthetic rate for a 
specified set of light conditions, nutrient concentration, 
and temperature is a product of the photosynthetic potential 
(P*RAD) and the maximum possible rate (NTLIM), The second 
term in the phytoplankton equation represents the loss of 
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biomass from respiration; 



RESP2 = RH020 ♦ X2 (29) 

RESP2 = flux of biomass expended _ 

by cell respiration (g C/m^/day) 

RH020 = rate of respiration per ^ 

unit of phytoplankton (day’' ) 

X2 = standing stock of phytOr- 2 

plankton (g C/xn), 

Phytoplankton respiration rate is defined as a function of 

temperature following Riley (,1946); 



RH020 = RH0Z2 * expCR*TEMP) (30) 

RH020 = rate of respiration per 

unit of phytoplankton (day’" ) 

RH0Z2 = rate of respiration at - 

0°C, (day’'-^) 

R = constant, rate of change 

in respiration rate with ^ 

temperature (°C" ) 

TEMP = mixed layer temperature (®C), 



The final expression in the phytoplankton equation represents 
the loss of biomass to grazing by herbivorous zooplankton: 



GRAZ = TAU23 * X3 (31) 

GRAZ = flux of biomass from 

phytoplankton to 2 

zooplankton (g C/m^/day) 

TAU23= zooplankton grazing rate (day~^) 

X3 = standing stock of 2 

zooplankton (g C/m ). 

A saturation grazing rate is defined by a Michaelisr-Menten 

type equation (Fig. 3) following Walsh and Bass (1971): 
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(32) 



TAU23 - MAXG*CX2MGrvX2MIN)/(KG+(X2MG-X2MIN)) 



X2MG 



KG 

3, 



zooplankton grazing rate 


(day"^) 


maximum specific grazing rate 


(day"^) 


concentration of phytoplankton 
in a mixed layer 


(mg C/m^) 


minimum concentration of 
phytoplankton below which 
grazing will not occur 


(mg C/m^) 


half saturation constant 


(mg C/m^) . 



Zooplankton Equation 
The zooplankton equation is; 
dX3 



dt 



= GRAZ - RESP3 r- VOID ^ LOSS 



which in words is: 



The change The input from The loss 

in zooplank- grazing on of biomass 

ton biomass (= )phytoplankton(") from respi- 
with time biomass ration 



The loss of 
biomass by 
excretion 



The loss of 
biomass by 
(-) predation and 
mortality 



The first term (GRAZ) in this equation has been defined 
(see eq, 31). The second term is a temperature dependent 
respiration term similar to that for phytoplankton: 



RESP3 = RH030 * X3 




(34) 


RESP3 = 


flux of biomass during 
respiration 


(g C/m^/day) 




RH030 = 


rate of respiration per 
unit of zooplankton 


(day“^) 




X3 


standing stock of 
zooplankton 


(g C/m^). 
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The equation for the respiration rate is: 



RH030 = RH0Z3 * exp(R*TEMP) (35) 

RH030 = rate of respiration per ^ 

unit of zooplankton (day'~ ) 

RH0Z3 = rate of respiration at 0°C. (day"^) 

R = rate of change in 

respiration rate as a - 

function of temperature (“O'" ) 

TEMP = mixed layer temperature (®C). 

The loss of biomass by excretion is: 



VOID = TAU30 * X3 



(36) 



2 

VOID = flux of voided biomass (g C/m /day) 

TAU30 = rate of excretion per 

unit of zooplankton ^ 

biomass (day*""^) 

X3 = standing stock of 2 

zooplankton (g C/m ), 

Based on data collected by Harvey et . al. (1935) 
and his own observations, Riley (1947) concluded that 
utilization (ratio of assimilation to ingestion) of grazed 
phytoplankton was very nearly unity when algal concentrations 
were low. During such periods, therefore, minimal carbon 
biomass would be voided. However, with the onset of the 
spring flowering, more food was available than could be 
used, so that quantities were voided in a semidigested 
condition. To approximate these observations, the fraction 
of biomass voided has been treated as a function of the 
grazing rate: 



TAU30 = (0.8*(TAU23/MAXG)) * TAU23 



(37) 




■"TA 



i 



TAU30 = rate of excretion per unit ^ 

of zooplankton biomass (day^"^) 

TAU23 = zooplankton grazing rate (day”^) 

MAXG = maximum specific grazing rate (day”^). 

When grazing rates are low (.025 day~^), the bulk of 
the ingested phytoplankton is assimilated (92%) and only 
8% is voided. During the bloom conditions observed in 
March 70% of the grazed biomass was assimilated with 30% 
being voided. These values are in general agreement with 
Conover’s (p. 113, in Parsons and Takahashi, 1973) conclusion 
that, in general, zooplankton assimilation efficiencies 
range from 60 to 95%. 

The final term in the zooplankton equation represents 
losses of zooplankton biomass by predation and natural 
mortality. 



LOSS = LAM30 * X3 (38) 

LOSS = flux of biomass due to 2 

predation and mortality (g C/m'^/day) 

LAM30 = rate of loss per unit of 

zooplankton (day” ) 

X3 = standing stock of 2 

zooplankton (g C/m^), 

The rate of loss has been estimated by this author as a 

constant 1% of the standing stock per day. 



LAM30 = L (39) 

L = fraction of zooplankton lost - 

per day (day” ). 
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TABLE I 



VALUES FOR EQUATION CONSTANTS 
AND THEIR RESPECTIVE SOURCES 



NAME 


DEFINITION 


VALUE 


REFERENCE 


B 


MAXN at 0°C 


0,50 


estimate 


C 


phosphate : carbon ratio 


0,774 


Riley (1965) 


DXl 


nutrient concentration 


3,0 


estimated from 
data 


EKL 


depth of Ekman layer 


50,0 


Smith (1968) 


FAVL 


nutrient regeneration 
coefficient 


0,8 


Butler , et al , , 
(1970) 


KG 


half saturation constant 
for zooplankton grazing 


80.0 


McAllister (1970) 


KN 


half saturation constant 
for phytoplankton grazing 


0.6 


Ichimura, (in 
Parsons and 
Takahashi, 1973) 


L 


zooplankton loss rate 


0.01 


estimate 


M 


coefficient of exchange 


0.12 


estimate 


MAXG 


maximum grazing rate 


0.25 


McAllister (1970) 


P 


photosynthetic constant 


2.5 


Riley (1946) 


R 


respiration coefficient 


0.069 


Riley (1946) 


RHOZ2 


phytoplankton respiration 
at 0°C. 


0,0175 


Riley (1946) 


RHOZ3 


zooplankton respiration 
at 0®C. 


0.019 


Riley (1947) 


S 


rate of change of MAXN 
with temperature 


0.13 


estimate 


XIMIN 


minimum nutrient uptake 
concentration 


0.0 


Dugdale (1967) 


X2MIN 


minimum phytoplankton 
concentration for grazing 


16.0 


McAllister (1970) 



42 



IV. PROCEDURE 



A. SIMULATION 

The system of differential equations was integrated 
numerically using the System/360 Continuous Simulation Model- 
ing Program (S/360 CSMP) available on the Naval Postgraduate 
School IBMr-360 computer. A brief discussion of S/360 CSMP 
is given in the Appendix. A more complete description is 
available in the IBM User'‘s Manual (1970). 

B, DATA COLLECTION 

Theoretical results, generated from observed environmental 
coiditions (i.e., radiation, temperature, etc,), were 
compared with observed seasonal variations in phosphate 
concentration and zooplankton biomass carbon. Data was 
obtained on a series of cruises conducted aboard the R/V 
ACANIA from February through November, 1974, Dissolved 
reactive phosphate was measured using a Technicon II Auto- 
analyser. Concentrations in the mixed layer were determined 
by computing a "weighted mean" of measured values within the 
layer. The "weight" for each value was assigned according to 
the interval between sample depths. 

Zooplankton samples were collected in a number three 
(333 vim mesh) net with a mouth diameter of c.a. 25 cm. The 
carbon content of each sample was determined using a high 
temperature dry combustion thermal conductivity method 
(Traganza, Radney, and Grahm, 1975). A complete listing 
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f 




of cruise data including phosphate, temperature, and carbon 
values is available, (Traganza, 1975). 

C. PARAMETER TESTING 

The response of the model to different environmental 
conditions was examined by generating simulations using 
hypothetical values for environmental functions. The values 
of incident radiation (RADI), mixed layer temperature (TEMP), 
mixed layer depth (Z) and upwelling velocity (W) were varied 
10% above and below their ’’standard value” (the observed or 
literature value). A seperate simulation was computed for 
each variation, leading to eight seperate cases. The 
variations tested are indicated in figure (4). The O's 
represent the standard or observed value, the +’s indicate 
a 10% increase, and the X's a 10% decrease. 
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V , RESULTS 



A. COMPARISON OF SIMULATION WITH OBSERVED DATA 

The initial simulation results are shown in Fig. (5), 
Theoretical cycles appear as asterisks and the observed 
cruise values are indicated by X's. The average difference 
between simulated and observed nutrient concentrations is 
66.4%. A 94% average error exists between the simulated 
and observed standing stock of zooplankton. Although the 
absolute error is large, the trends of the data and simulation 
are quite similar. The principle source of error in the 
phosphate cycle was a 20 day variation between the simulated 
and observed nutrient minimum. The zooplankton simulation 
error was caused by the failure of the model to match the 
apparently low biomass values observed in late March and 
early September. 

Seasonal changes in mean phosphate concentration in the 
mixed layer were characterized by a steady increase during 
winter and spring months from an initial value of 1.5 yg-atom/l 
to a maximum 1.86 yg-atom/1. The simulated concentration 
was assigned an initial condition of 1.5 yg-atom/1 and subse- 
quently increased to 1.99 yg-atom/1. In June a sharp decline 
in concentration occurred both in the data and in the simula- 
tion. A minimum value of 0.33 yg-atom/1 was observed in early 
August. The simulated minimum (0.61 yg-atom/1) occurred in 
late August. Both the theoretical and the measured 
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concentrations increased steadily from August through 
December culminating at 1,0 pgr^a tom/1. 

The differences between observed and simulated maximum 
and minimum nutrient concentrations indicated that the 
function controlling nutrient uptake, NTLIM, (eq. 27) did 
not fit existing conditions. Testing various values of the 
constants KN, C^q . 27) S and B (eq. 28) indicated that a 
reduction in the value of KN from 0,6 pg^atom/1 to 0,3 
yg-atom/1 produced a more accurate correspondence between 
measured and theoretical maximum and minimum phosphate 
concentrations (Fig. 6). 

Phytoplankton data was not available for comparison with 
theoretical values. However, the sharp decline in nutrient 
concentrations during June and July suggests the presence 
of substantial phytoplankton stocks during this period. The 
simulated bloom, which peaks in July, correlates with the 
observed decline in phosphate. The phytoplankton maximum 
does not, however, precede the initial zooplankton peak as 
is normally observed. 

The cycle of zooplankton biomass was characterized by 

2 

two distinct maxima. The first peak (1.05 gC/m ) was 

2 

observed in late July. A simulated maximum (1.02 gC/m ) 

occurred during the same period. This initial peak was 

followed by a sharp decline of measured biomass. A minimum 
2 

of 0.18 gC/m was observed in early September. A lesser 

2 

reduction occurred in the theoretical values (0.78 gC/m ). A 
rapid increase in the standing stock followed the September 
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2 

minimum. An overall maximum of 1,85 gC/m was observed in 

mid-December. Simulation values also increased producing 

2 

a lesser maximum of 1,58 gC/m , 

The simulated energetic balance observed in the zooplankton 
varied with changes in phytoplankton concentration and 
temperature, During the onset of the phytoplankton bloom 
in mid^March, zooplankton assimilated 70% of the ingested 
carbon biomass and voided 30%, Of that carbon which was 
assimilated 60% was lost through respiration. In mid-July, 
when phytoplankton stocks had peaked, the high algal 
concentration resulted in near maximum grazing rates. Because 
more biomass was being ingested than could be used, 77% was 
voided and 23% assimilated. Due to higher mixed layer 
temperatures, 75% of the 23% assimilated biomass was required 
for respiration. 

Calculation of the annual average of monthly zooplankton 
gains and losses indicated 65% of the ingested phytoplankton 
biomass was assimilated, 35% was voided. Of the 65% 
assimilated biomass, 86% was respired leaving 14% (9% of 
that assimilated) for growth. Of this remaining biomass, 

75% was lost to predation and natural mortality. 

B. EFFECTS OF VARIATIONS IN ENVIRONMENTAL PARAMETERS 

Additional simulations were generated to determine the 
response of state variables to changes in the environmental 
parameters. Table II lists the percent change each modified 
parameter produced in the phytoplankton and zooplankton 
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TABLE II 



CHANGE IN MAXIMUM BIOMASS OBSERVED 
WITH +10% CHANGE IN FORCING FUNCTIONS 

PHYTOPLANKTON ZOOPLANKTON 

RADI, r.10% 

RADI, +10% 

LAYER, -10% 

LAYER, +10% 

UPWEL, -10% 

UPWEL, +10% 

TEMP, -10% 

TEMP, +10% 



-13% 


-07% 


+ 11% 


+10% 


+01% 


+02% 


^02% 


▼-03% 


-01% 


-01% 


00% 


00% 


-10% 


+85% 


+06% 


-65% 
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standing stock maxima, Complete simulation results are 
given in Figures (7 - 10). The 10% increase in incident 
radiation produced a nearly equivalent gain in the maximum 
biomass of both phytoplankton and zooplankton (Fig. 7). 
Increased radiation produced higher rates of photosynthesis 
and larger stocks of phytoplankton. These phytoplankton 
stocks in turn supported increased grazing by zooplankton. 

Reducing the layer depth produced a slight increase in 
biomass. Confining the plants in shallower water apparently 
exposed the phytoplankton stock to higher average radiation 
intensities and increased algal availability to zooplankton. 

Surprisingly, a 10% change in upwelling velocity had 
negligible effects on plankton biomass. Increasing velocities 
did produce higher nutrient concentrations (Fig. 9). However, 
this change in concentration did not yield a similar increase 
in production. This result stems from the non-linear relation 
ship between nutrient concentration and photosynthetic rate 
(Fig. 2). Above 1.0 yg-atom/1, changes in nutrient concen- 
tration have a small effect on photosynthesis. Since the 
variations in the simulated nutrient concentrations were 
observed principally at concentrations above 1.0 yg-atom/1, 
these fluctuations had a minimal effect on plankton production 
The model's response here agrees with Odum's (1971) contention 
that nutrients are not a limiting factor for phytoplankton 
production in regions of upwelling, 

The largest variation in zooplankton biomass (-65 to +85%) 
was the result of a +10% variation in the mixed layer 
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temperature. Temperature changes directly affect the 
respiration rate of the zooplankton. Increased temperatures 
produced higher respiration rates and a net loss of biomass. 
Conversely, lower temperatures produced lower respiration 
rates and greater standing stocks. Although the change in 
daily respiration losses is only slightly affected by temper- 
ature (7%/°C @ 15®C), the cumulative affect over a period of 
months appears significant Csee Fig, 10), 
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VI. DISCUSSION 



The primary environmental factor influencing simulated 
phytoplankton stocks was incident radiation, followed 
closely by temperature effects (see Table II), The greatest 
fluctuations in zooplankton biomass resulted from changes in 
the temperature regime. The response of zooplankton stocks 
to different thermal conditions suggests a hypothesis 
concerning the rapid September (warm water period) decrease 
in zooplankton stocks. 

An examination of mixed layer temperatures and predicted 
upwelling velocities for July and August (Fig, 4) indicate 
that this period was characterized by a reduction in upwelling 
(less cold deep water moving towards the surface) which 
resulted in higher mixed layer temperatures. These higher 
temperatures may have increased zooplankton respiration 
sufficiently to produce a net decrease in biomass. In late 
fall, when lower temperatures returned with algal stocks 
remaining sufficiently large to permit grazing, a rapid 
zooplankton increase occurred. 

Certainly, other processes contribute to changes in zoo- 
plankton stocks. The late summer decline may have been the 
result of intense predation during this period. Since the 
model assumes a constant loss to predation and natural 
mortality, variations in predation have not been considered. 
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A second possibility is that water masses containing 
different size plankton stocks may have appeared in Monterey 
Bay in late summer and winter. The observed changes in the 
bay have traditionally been described in terms of three 
oceanographic periods (Skogsberg, 1936; cited in Sverdrup 
et al., 1942), The "Davidson Current period" (relatively 
warm water) is normally observed in winter, the "upwelling 
period" (cold water phase) from early spring through summer, 
and the "oceanic period" (warm water phase), from late summer 
through fall, A comparison of these dates with the observed 
zooplankton fluctuations indicates that the first maximum 
might be associated with the upwelling period, the minimum 
with the oceanic, and the second maximum with the Davidson 
Current period. 

Although changing current patterns clearly exist in 
Monterey Bay, the role of advection in controlling plankton 
populations may be neglected if (for the time and space scale 
of interest) the region can be assumed horizontally homogeneous. 
The total rate of change of a variable "S" as given by the 
material derivative (^) is: 



= IS ^ ^ 

Dt 9t ^ 9x ^ 9y 



( 1 ) 



.DS. 



The total material derivative (p^) will equal a simple time 
dS 

derivative (^) if the horizontal velocities (u, v) are zero 
or if the region is horizontally homogeneous, i.e., the 
spatial' derivatives -|^) art? zero. 
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The argument has been made by Bolin and Abbott (1962) 
that the hydrographic and biological features observed in 
Monterey Bay are characteristic of a much wider geographic 
region. These investigators cited evidence that upwelling 
in the bay is paralleled by upwelling along most of the West 
Coast of the United States north of Point Conception. Their 
temperature observations also corroborated closely with 
the more extensive data reported in CalCOFI Report, VII, 1960. 
In addition similar phytoplankton observations were made at 
a station as far distant as La Jolla, Their conclusion was 
that the assumption of horizontal homogeneity is appropriate 
for this region, 
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VII, CONCLUSIONS 



The ecosystem model simulated with reasonable accuracy 
the observed seasonal changes in phosphate concentration 
and zooplankton biomass. The model inputs used to generate 
the theoretical results consisted of parameters defining those 
characteristics of the physical environment which effect bio- 
logical productivity, e.g., solar radiation, temperature, 
layer depth, and nutrient input. 

Analysis of the model's sensitivity to changes in these 
environmental parameters identified the possible importance 
of thermal conditions in regulating the growth of zooplankton 
populations. These theoretical results suggest a hypothesis 
related to upwelling. The importance of upwelling as a 
source of nutrients has long been acknowledged. Simulation 
results indicate that upwelling may also enhance the growth 
of zooplankton stocks by bringing colder deep water to the 
surface thereby reducing zooplankton respiration requirements. 

Of course, the model can only suggest the possible signif- 
icance of this cooling effect. Actual verification of the 
hypothesis will require a precise knowledge of zooplankton 
metabolism, as well as additional measurements of phytoplankton 
and zooplankton biomass. A more detailed simulation study 
of the effects of upwelling might best be made using a two- 
dimensional model to analyze variables as a function of 
depth as well as time. 
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FIGURE 1. Ecosystem Compartment Model and Flux Equations. 
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FIGURE 2. Function NTLIM (maximum photosynthetic rate). 
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FIGURE 3. Function TAU23 (zooplankton grazing rate). 
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FIGURE 6. Simulation Results (half-saturation constant, 
KN = 0.3). 
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FIGURE 7, Simulation Results (incident radiation, RADI, +10%) 
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FIGURE 8. Simulation Results (mixed layer depth, Z, +10%) 
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FIGURE 9, Simulation Results (upwelling velocity, W, +10%), 
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APPENDIX A 



A BRIEF DESCRIPTION OF THE S/360 CONTIUOUS SIMULATION 
MODELING PROGRAM (CSMP) 

The S/360 CSMP program is devided into three segments, 
identified by the statement labels INITIAL, DYNAMIC, and 
TERMINAL. The INITIAL segment is intended for the specifi- 
cation of initial conditions and for defining constants and 
parameters. The DYNAMIC segment includes the complete 
description of the system, together with any other computa- 
tions desired during the run. The TERMINAL segment is used 
for those computations desired after the completion of each 
run. The INITIAL and TERMINAL segments are optional, but 
the DYNAMIC segment is mandatory, 

Structure Statements 

Structure statements define the model to be simulated. 
In general, rules for structure statements follow those used 
in FORTRAN. A library of special CSMP functions is also 
available. One such function used in the present model is 
the INTGRL function. This function computes numberically the 
integral of a function ”X." It is expressed in the general 
form : 

Y = INTGRL (IC, X) 
where the initial condition is; 
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IC = Y(P), 



This CSMP statement is equivalent to | 

Y = J X dt + IC 

O 

Data Statements 

Data statements are used to assign numerical values 
to constants, parameters, and initial conditions. These 
data cards are specified by the labels CONSTANT, PARAMETER, 
and INCON respectively. If a specific parameter is to be 
tested over a range of values a sequence of simulation runs 
may be designated by enclosing several values of the variable 
in parentheses. For example; 

PARAMETER X = (5.0, 5.1, 5.2), 

specifies three simulation runs. The value of X will be 5.0 
in the first run, 5.1 in the second, and 5.2 in the third. 

A linear function connecting data points may be 
specified by the function generator AFGEN. The set of data 
points is identified using the FUNCTION label: 

FUNCTION UPWELL = (O. , .10), (10., .50). 

The points are given in sequence, X^-coordinate (time) first 
followed by the y-coordinate . The program connects the data 
points with a linear function when the AFGEN function is 
called; 



W = AFGEN(UPWELL, TIME). 
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Control Statements 



Certain operations related to translation, execution, 
and output segments of the program are specified by special 
control statements; 

These three labels identify the major segments 
of the CSMP program. 

This statement marks the completion of the model's 
structural description. 

These labels determine whether a sequence of cards 
is to be machine sorted into a correct sequence 
or executed in the order given. 

This card follows the last END statement in the 
program. 

This statement denotes the end of a job and must 
follow the stop card. 

This label is used with the following CSMP 
specified variables. 

Output print increment. 

Print-plot output print increment. 

The maximum time for simulation. 

This is only a partial listing of the functions and 
labels used in the S/360 CSMP language. For a complete 
description see the IBM User's Manual (1970). 



INITIAL 

DYNAMIC 

TERMINAL 

END 

SORT 

NOSORT 

STOP 

ENDJOB 

TIMER 

PRDEL 

OUTDEL 

FINTIM 
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FUNCTION SOLAR=(OO.t.l23), (46.r.l58), (74. ,.224), 
(105., .286), ( 135. ,.324), ( 166. ,.335), (196. ,.31 1 ) 
( 227., .278), (258. ,.235), ( 288. ,. 181), (319. ,.141) 
(349. ,.11 7), ( 365. ,. 119) 
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